Model independent study of the Dirac structure of the nucleon-nucleon interaction 



O. Plohl, C. Fuchs, and E. N. E. van Dalen 

Institut fur Theoretische Physik, Universitat Tubingen, 
Auf der Morgenstelle 14, D-72076 Tubingen, Germany 

Relativistic and non-relativistic modern nucleon-nucleon potentials are mapped on a relativistic 
operator basis using projection techniques. This allows to compare the various potentials at the 
level of covariant amplitudes were a remarkable agreement is found. In nuclear matter large scalar 
and vector mean fields of several hundred MeV magnitude are generated at tree level. This is found 
to be a model independent feature of the nucleon-nucleon interaction. 
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A fundamental question in nuclear physics is the role 
which relativity plays in nuclear systems. The ratio of 
the Fermi momentum over the nucleon mass is about 
k^/M ~ 0.25 and nucleons move thus with maximally 
about 1 /4 of the velocity of light which implies only mod- 
erate corrections from relativistic kinematics. However, 
there exists a fundamental difference between relativis- 
tic and non-relativistic dynamics: a genuine feature of 
relativistic nuclear dynamics is the appearance of large 
scalar and vector mean fields, each of the magnitude of 
several hundred MeV. The scalar field E5 is attractive 
and the vector field S M is repulsive. In relativistic mean 
field (RMF) theory, both, sign and size of the fields are 
enforced by the nuclear saturation mechanism [lj. 

At nuclear saturation density po — 0.16 fm -3 the em- 
pirical fields deduced from RMF fits to finite nuclei are 
of the order of S s ~ -350 MeV and E Q ~ +300 MeV 
21 (E is the time-like component of £ M ). The sin- 
gle particle potential in which the nucleons move orig- 
inates from the cancellation of the two contributions 
U s . p , ~ E + E s ~ -50 MeV which makes it difficult 
to observe relativistic effects in nuclear systems. There 
exist, however, several features in nuclear structure which 
can naturally be explained within Dirac phenomenology 
while models based on non-relativistic dynamics have dif- 
ficulties. Best established is the large spin-orbit splitting 
in finite nuclei. Also the so-called pseudo-spin symmetry, 
observed more than thirty years ago in single-particle lev- 
els of spherical nuclei, can naturally be understood within 
RMF theory as a consequence of the coupling to the lower 
components of the Dirac equation Q . 

A connection to Quantum-Chromo-Dynamics (QCD) 
as the fundamental theory of strong interactions is estab- 
lished by QCD sum rules [j]. The change of the chi- 
ral condensates (qq), (q^q) in matter leads to attractive 
scalar and repulsive vector self-energies which are aston- 
ishingly close to the empirical values derived from RMF 
fits to the nuclear chart. Also relativistic many-body cal- 
culations 0, 0, yield scalar/vector fields of the same 
sign and magnitude as obtained from RMF theory or, 
alternatively, from QCD sum rules. Moreover, Dirac- 
Brueckner-Hartree-Fock (DBHF) calculations p| agree 



even quantitatively surprisingly well with the QCD moti- 
vated approach of Ref. |9| where chiral fluctuations from 
pion- nucleon dynamics were considered on top of the chi- 
ral condensates. 

These facts suggest that preconditions for the existence 
of large fields in matter or, alternatively, the density de- 
pendence of the QCD condensates, must already be in- 
herent in the vacuum nucleon-nucleon (NN) interaction. 

Relativistic one-boson-exchange potentials (OBEP), 
e.g. Bonn (A,B,C), ^| provide fits with fair preci- 
sion to NN scattering data. High precision fits (with 
X 2 /datum ~ 1.01) are provided by QBE type potentials 
such as the relativistic CD-Bonn [ll| or the Nijmegen po- 
tentials ^3] where fits are performed separately in each 
partial wave. In contrast to e.g. the spectator model 
of 13] the 'standard' OBE potentials are based on the 
no-sea approximation which excludes explicit excitations 
of anti-nucleons but absorbs such contributions into the 
model parameters, in particular into a large uj coupling. 
The present investigations are restricted to the standard- 
type models. 

The long-range part of the interaction is generally me- 
diated by the one-pion-exchange (OPE) while the scalar 
intermediate range attraction is mainly due to correlated 
two-pion-exchange. The short-range part, i.e. the hard 
core, is dominated by the heavy vector meson exchange, 
i.e. the isoscalar u> and the isovector p meson. Widely 
used in nuclear physics are, however, also high precision 
non-relativistic empirical potentials such as the Argonne 
potential AVis 14] • A systematic connection to QCD is 
established by chiral effective field theory (EFT). Up to 
now the two-nucleon system has been considered at next- 
to-next-to-next-to-leading order (N 3 LO) in chiral pertur- 
bation theory 0,0]- I n sucn approaches the NN poten- 
tial consists of one-, two- and three-pion exchanges and 
regularizing contact interactions which account for the 
short range correlations. The advantage of EFT com- 
pared to the more phenomenological OBE potentials is 
the systematic expansion in terms of chiral power count- 
ing. 

A better understanding of the common features of the 
various approaches is essential in order to arrive at a 
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more model independent understanding of the nucleon- 
micleon interaction, in particular since all well estab- 
lished interactions fit NN scattering data with approx- 
imately the same precision. Substantial progress was re- 
cently achieved by the construction of a universal low en- 
ergy NN potential Vl ow k based on renormalization group 
techniques ^t}. However, like the EFT potentials V\ ovl k 
is not covariantly formulated and has therefore not been 
used in relativistic nuclear structure calculations. 

The present work applies projection techniques to map 
the various potentials on the operator basis of relativis- 
tic field theory which is given by the Clifford algebra 
in Dirac space. This allows to identify the different 
Lorentz components of the interaction and to calcu- 
late the relativistic self-energy operator in matter. The 
philosophy behind this approach is based on the fact, 
that any NN interaction, whether relativistic or not, is 
based on a common spin-isospin operator structure which 
invokes certain scales: long-range spin-isospin depen- 
dent forces, essentially given by the one-pion exchange, 
short- and intermediate-range spin-independent interac- 
tions, short-range isoscalar spin-orbit interactions and 
quadratic isovector spin-orbit interactions. 

In a covariant formulation the spin-isospin structure 
of the interaction is constrained by the symmetries of 
the Lorentz group. E.g. the Born scattering matrix of 
covariant OBE potentials is given by the sum over the 
corresponding scalar, pseudoscalar and vector mesons a 

V(q',q) = ^'9)4 2) B a («/-9)rW , (1) 

a—s^ps,v 

where q^ and q'^ are the cm. momenta of the incom- 
ing and outgoing nucleons. The meson propagator is 
denoted by D a , T a are form factors applied to the ver- 
tices T Q . Isospin factors (t\-t-i) are suppressed in JQ. 
The Dirac structure of the potential is contained in the 
meson-nucleon vertices 

r s = g s l, T ps = g ps — ?7 5 , i\, = g v ^ + — .(2) 

For the pseudoscalar mesons 7r and r\ a pseudovector cou- 
pling is used and the lu meson has no tensor coupling, i.e. 

= 0. The potential V(q!,q), i.e. the OBE Feynman 
amplitudes are obtained by sandwiching V between the 
incoming and outgoing Dirac spinors. 

From a low energy expansion of the OBE amplitudes 
one obtains the representation of non-covariant poten- 
tials 

nq',q) = E[^ + y - T i- T 2] o a (3) 

where the most important operators (assuming identicle 
particles) are 

0\ = 1, 02= erx-o-2, 

3 = ( < T 1 -k)( < T 2 -k)-ik 2 ( t Ti- CT2 ), 

4 = |(o-i + <T 2 )-n, 5 = (<ri-n)(o- 2 -n), (4) 



with k = q' — q, P = i(q' + q) and n = P x k. 

The modern non-covariant Nijm I and II potentials 
are constructed from approximate OBE amplitudes and 
thus based on the operator structure (@J . However, since 
they are separately fitted in each partial wave they are 
often referred to as phenomcnological models. The Ar- 
gonne potential AVis fl^t is purely phenomenological in 
the sense that the OBE picture is released. Based on the 
symmetries of (@J , the intermediate and short range part 
is parameterised by phenomcnological functions V a and 
only one-pion-exchange is contained explicitly. The EFT 
potentials, containing one- and multi-pion exchange ex- 
plicitely, are even more rigorous since most part of the 
nuclear repulsion is carried by regularizing counter terms. 
Such an approach is justified by the fact that heavy me- 
son (p, uj) exchange cannot be resolved up to the mo- 
mentum scale of about 400 MeV which is constrained by 
NN scattering data. We apply the EFT Idaho potential 
[l^j l which fits NN scattering data with similar quality as 
AVis or the Nijmegen potentials. The same philosophy 
is behind V\ ovr k which can be viewed as the condensa- 
tion of the various formulations to a model independent 
result. 

Although based on the same operator structure as 
OBEPs, it is not straightforward to fix the Lorentz char- 
acter of the various pieces of the non-covariant interac- 
tions AVis, EFT (Idaho) and Mow k- 

However, any two-body amplitude can be represented 
covariantly by Dirac operators and Lorentz invariant 
amplitudes [l8j. A relativistic treatment invokes au- 
tomatically the excitation of anti-nucleons. However, 
NN scattering, in both, relativistic and non-relativistic 
approaches is restricted to the positive energy sector 
and neglects the coupling to anti-nucleons. As a con- 
sequence one has to work in a subspace of the full 
Dirac space where on-shell two-body matrix elements 
can be expanded into five Lorentz invariants. A possi- 
ble choice of a set of five linearly independent covariant 
operators are the scalar, vector, tensor, axial-vector and 
pseudo-scalar Fermi covariants r m = {S, V, T, P, A} with 
S = 101, V = 7^®7 M , T = <7^<g>o>, P = 75^75, A = 
757^0757^. Working with physical, i.e. antisymmetrized 
matrix elements V A , one has to keep in mind that direct 
and exchange covariants T m (where Dirac indices of par- 
ticles 1 and 2 are interchanged) are coupled by a Fierz 
transformation. Since the low momentum part of the 
NN interaction is totally dominated by the pseudovector 
one-pion exchange (OPE) it is preferable to choose an 
operator basis where the pseudovector exchange is com- 
pletely separated from the remaining operator structure. 
This can be achieved by the following set of covariants 
originally proposed by Tjon and Wallace 0] 

r m = {S,-S,(A-A),PV,-PV} . (5) 

PV and PV are the direct and exchange pseudovector 
covariants, analogous to the pseudoscalar covariant P, 
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however, with 75 replaced by — $)/2M^. Thus the 
on-shell (|q| = |q'|) scattering matrix is given by 



^(q',q) = £s™(|q|,< 



(6) 



where 9 is the cm. scattering angle and I = 0, 1 the 
isospin channel. For the Hartree-Fock self-energy it is 
sufficient to consider 9 — when antisymmetrized ma- 
trix elements are used since 9 = n contains then only 
redundant information. The transformation of the Born 
amplitudes from an angular-momentum basis onto the 
covariant basis JHJ is now standard and runs over the 
following steps: \LSJ) — * partial wave helicity states — > 
plane wave helicity states — > covariant basis. The first 
two transformation can e.g. be found in Refs. |H^. The 
last step has to be performed numerically by matrix in- 
version FlEoll. 
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FIG. 1: Isospin-averaged Lorentz invariant amplitudes 
<7m(|q|,0) for the different NN potentials after projection on 
the Dirac operator structure. The pseudovector representa- 
tion of the relativistic operator basis is used. 

Fig. Hshows the resulting isospin-averaged amplitudes 
9m = k(giT° + Sff^T 1 ) at 6 = for the various potentials 
up to 400 MeV. The four amplitudes gs , <?§, 5a and 
g^y are very close for the OBEPs Bonn A, CD-Bonn, 
Nijmegen I and II, and the non-relativistic AVis, Idaho 
and Vjow k- The direct pseudovector amplitude gpy falls 
somewhat out of systematics. This amplitude is, how- 
ever, of minor importance since it does not contribute 
to the Hartree-Fock self-energy (J7J and to the single 
particle potential. The dominance of the OPE at low 
q| is reflected in the pseudovector exchange amplitude 
$pv wmc h is at small |q| almost two orders of magni- 
tude larger than the other amplitudes. For the covariant 
OBEPs Bonn and CD-Bonn and the OBE-like Nijmegen 



potentials a coincicdence at the level of covariant am- 
plitudes is not completely unexpected. However, it is re- 
markable that this agreement transfers to the EFT Idaho 
potential and also to V\ ow k which are of completely dif- 
ferent character and theoretical background. 

With the covariant amplitudes at hand, one is now able 
to determine the relativistic mean field in nuclear matter 
by calculating the relativistic self-energy E in Hartree- 
Fock approximation at tree level. We are thereby not 
aiming for a realistic description of nuclear matter sat- 
uration properties which would require a self-consistent 
scheme. Moreover, short-range correlations require to 
base such calculations on the in-medium T-matrix rather 
than the bare potential V [2]]]. This leads to the rela- 
tivistic Dirac-Brueckner-Hartree-Fock scheme which has 
been proven to describe nuclear saturation with quanti- 
tatively satisfying accuracy 0, Q> HI 

The self-energy for the nucleon k follows from the inter- 
action matrix V by integrating over the occupied states 
q in the Fermi sea 



H a p(k,k F ) 



G»M V A (k,q) a ^ T . (7) 



The Dirac propagator G D {q) = [jf + M)2iti5{q 2 



M 2 )Q(q )Q(kF — |q|) describes the on-shell propaga- 
tion of a nucleon with momentum q and energy Eq — 

In isospin saturated 
a 



i/q 2 + M 2 inside the Fermi sea. 
nuclear matter the self-energy consists of a scalar E s 
time-like vector E Q and a spatial vector part E v 



E(fc, k F ) = E 8 (fe, k F ) - 70 E D (fc, k F ) + 7 ■ kE v (fc, fc F ), (8) 
which are given by 7] 
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Fig. |2 shows the tree level self-energy components 
obtained with the various potentials in nuclear matter 
at saturation density with a Fermi momentum k F = 
1.35 fra~ . All potentials yield scalar/vector fields E s 
and E Q of comparable strength: a large and attractive 
scalar field E s ~ — (450 -j- 400) MeV and a repulsive vec- 
tor field of E ~ +(350 400) MeV. The spatial compo- 
nents kE v agree as well. As known from self-consistent 
DBHF calculations 0, , the spatial vector self-energy 
is a moderate correction to the large scalar and time-like 
vector components E s and E Q . This is already the case at 
tree level. The right part of Fig. shows the single par- 
ticle potential U s . p .(k, k F ) = ME S /E k - E Q + E v k 2 /#k, 
defined as the expectation value of the self-energy. It 
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FIG. 2: Tree level results for nuclear matter at &f = 1.35 fm~ for the various potentials: scalar and vector self-energy 
components (left), spatial vector self-energy component (middle) and single particle potential Us.p. (right). Stars in the right 
panel denote a non-relativistic calculation of U s . p .. 



reflects the well known fact that various two-body po- 
tentials are rather different although they are phase-shift 
equivalent, i.e. they describe NN scattering data with 
about the same accuracy when iterated in the Lippmann- 
Schwinger equation |2l(. The differences at tree level 
are mainly due to differences in the short-range part of 
the interaction 01 ■ Fig- 121 includes also the results 
from a 'non-relativistic' calculation of U a . p . where the 
partial wave amplitudes are directly summed up. The 
non-relativistic and the relativistic calculations show an 
excellent agreement which demonstrates the accuracy of 
the applied projection techniques. 

We presented a model independent study of the Dirac 
structure of the nucleon-nucleon interaction. The poten- 
tials were projected on a relativistic operator basis in 
Dirac space using standard projection techniques which 
transform from a partial wave basis, where the potentials 
are originally given, to the basis of covariant amplitudes. 
The different approaches can now be compared at the 
level of these covariant amplitudes and, moreover, this 
allows to calculate the relativistic self-energy operator 
in nuclear matter. For both, the covariant amplitudes 
and the tree-level Hartree-Fock self-energy, we observe a 
remarkable agreement between the OBEPs (Bonn, CD- 
Bonn, Nijmegen), the phenomenological AVis potential 
and the EFT based Idaho and Mow k potentials. The 
structure of the interaction enforces the existence of large 
scalar and vector fields as a model independent fact. The 
scale of these fields is set at the tree level. Although es- 
sential for nuclear binding and saturation, higher order 
correlations, i.e. Brueckner ladder correlations, change 
the size of the fields by less than 20% Q. The magnitude 
of the fields is similar to that predicted by relativistic 
mean field phenomenology, relativistic many body corre- 



lations and also by QCD sum rules. We conclude that 
the appearance of large scalar and vector fields in mat- 
ter is a general and model independent consequence of 
the vacuum structure of the NN interaction. Relativistic 
dynamics is therefore essential for nuclear systems. 

We like to thank D. Entem for providing the Idaho 
N3LO program package and H. Muther for providing his 
Hartree-Fock code. 
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